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• and their time discretization which preserves the constant area speed property. We 

■ establish various polygonal analogues of geometric formulas for moving boundaries 

■ and make use of the geometric formulas for our numerical scheme and its analysis 
' of general constant area speed motion of polygons. Accuracy and efficiency of our 
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1 Introduction 



Polygonal analogues of several moving boundary problems and their time discretization 



^ ■ are investigated in this paper. The polygonal motion is restricted within an equivalent 

Tj- ■ class of polygons. We introduce notion of polygonal curvature, which is consistent to 

•/^ ' polygonal analogues of geometric variational formulas. 

(yQ \ We propose a formulation of general area-preserving motion of polygonal curves by 



using a system of ODEs. The moving polygon belongs to a prescribed class of polygons, 
which is similar to the admissible class in the theory of crystalline motion by curvature. 
There are many articles about the crystalline curvature flow and asymptotic behavior of 
solutions [H El El El El Uni HH [IS, US, [E] , etc., which started from the pioneer works [2] and 
[TT] . Actually, if the initial curve is a convex polygon in a crystalline admissible class, then 
our polygonal curvature flow is nothing but the crystalline curvature flow. However, we 
consider more general polygonal moving boundary problems in wider admissible classes 
of polygons. 
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Based on the formulation of general polygonal moving boundary problems, we propose 
an implicit time discretization scheme with effective iteration scheme for the nonlinear 
system in each time step. It has second order accuracy and preserves the constant area 
speed property. In a fixed admissible class of polygons, we prove a convergence theorem 
of second order for our numerical scheme. 

On the other hand, it is expected that our polygonal analogue becomes a natural 
approximate solution of a smooth moving boundary problem if the number of edges is 
enough large. It is important and interesting application of our polygonal motion, but 
we do not touch on this issue in this paper. We only mention here that the crystalline 
algorithm and their convergence theorems for the motion by curvature are found in [5l El 
m El [m [121 [13] etc. 

The organization of this paper is as follows. Fundamental notation and formulas for 
polygons and polygonal motions are introduced in Section [2l In Section [3l a general 
initial value problem of polygonal motion in an equivalent class is considered, and the 
constant area speed condition is given. Several basic examples of polygonal motions such 
as the polygonal curvature flow and the polygonal advected flow are also presented. In 
Section |H an implicit scheme of Crank-Nicolson type and an iteration scheme for the 
nonlinear system in each time step are proposed. The proposed scheme inherits the 
constant area speed property and its second order convergence is proved in Theorem 14.51 
In Section [5l the accuracy of our numerical scheme is checked through various numerical 
simulations in comparison with the first order explicit Euler scheme. These simulation 
show that the second order scheme preserves the constant area speed property with high 
accuracy. 

2 Polygons and polygonal motions 

We give basic definitions and notation for the polygonal motion in an equivalent class of 
polygons. In particular, the polygonal curvature is introduced as a generalization of the 
crystalline curvature. We also collect their fundamental formulas and properties in this 
section. 

2.1 Polygons 

We define a set of polygons in M^: 

V := {T; r is a polygonal Jordan curve in R^}. 

In this paper, we assume that any two dimensional vector a; G is represented by a 
column vector, and we denote its transposed row vector by x"^. For F G P, the bounded 
interior polygonal domain surrounded by F is denoted by Q. For simplicity, we consider 
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Figure 1: Some quantities defined on Tj. 

the case that Q is simply connected, but many of the following arguments are valid in 
other geometrical situations. 

Let r G P be an A^-polygon. The N vertices of F are denoted by Wj G for 
j = 1,2,...,N counterclockwise, where Wq = and ioat+i = Wi. Hereafter we use 
the periodic numbering convention Fq = Fa? and Fat+i = Fi for any quantities defined on 
A^-polygon. 

For j = 1, 2, . . . , A^, the j-th edge between Wj_i and Wj is defined by 

r^. = {(1 - 0)wj^i + Owj- < < 1}, 

and its length is denoted by |Fj| := \wj — Wj^i\. The characteristic function Xj ^ L°^{V) 
for Tj is defined as 

xA^)--=\^. ""ll',^ (j = l,2,...,iV). 
10, a; G r \ Fj 

The outward unit normal on F^ is denoted by n^, and the outer angle at the vertex 
Wj is denoted by ipj G (— vr, vr) \ {0}. They satisfy cosipj = rij^i ■ rij. We also define the 
height of Tj from the origin by hj := Wj ■ rij = Wj^i ■ rij (see Figured]). Then the straight 
line including Tj is expressed by the equation rij ■ x = hj, and the vertices of F G P are 
given from {hj}j as 

Proposition 2.1 Under the above conditions, we have 

\Tj \ = aj„ihj_i + hjhj + ajhj^i (j = 1, 2, . . . , A^), (2.2) 
where aj := cosec ipj and hj := — cot f^j-i — cot ipj. 
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Proof. We define a unit tangent vector of Tj by tj := [wj — Wj_i)/\rj\. We remark tliat 
nj+i ■ tj = —rij ■ tj^i = siinpj and tj ■ tj^i = cosipj. Tlien, from (12. ip and tlie equality: 



nj \ / \ ( —rii ■ t 



we liave 



Since 

ir^l = tj-{wj-Wj^i) 

= ttjhj^i — {aj cos(pj + aj_i cos(y9j_i} hj + aj^ihj^i 
= (ijhj^i -\- bjhj -\- Qjj^ihj^ij 

we obtain the formula (12.21) . I 
The total length of F is given by 

N N N 

|r| := 5^ irjl = + bj + aj_i)hj = ^ rj^hj, (2.3) 

j=i j=i j=i 

where rjj := aj + bj + aj_i = tan((y9j/2) + tan(y9j_i/2). The area of interior domain Q is 
denoted by which is given by 

1 ^ 

l^l = 2E|r.l^r (2-4) 

i=i 

The above symbols are also written as nj = '^'^(r), aj = cij(r) and hj = hj{r) etc., 
provided we need to distinguish from quantities of the other polygons. 

2.2 Equivalence classes of polygons 

For two polygons T and S G P, we define an equivalence relation F ~ S. We say F ~ S, 
if their numbers of edges are same (let it be A^) and ^^(F) = ^^(S) for all j = 1, 2, . . . , 
after choosing suitable counterclockwise numbering for F and S. The equivalence class of 
F e P is denoted by V[T] := {S G P; S ~ F}. 

We fix an A^-polygon F* G P and let V* := V[T*]. For F and E in P*, we define the 
distance between them by 

d{T, S) := ^.^max^ \hj{T) - /ij(S)|. 
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Then, it is clear that {V*,d) becomes a metric space, since it is isometrically embedded 
in equipped with maximum norm | ■ |oo by the height function h defined on V*: 

h{T) := ih{T), h^iV), h^iV)) G {T E V*). 

We assume that vectors in are represented by row vectors. It is obvious that the 
image of the height function h(V*) is open in M^. 

For r G P* and e > 0, an e-ball in V* = V[T] with center F is denoted by 

B{T,e) := {S G V[T]; d(S,F) < e}. 
For an open set O G V* and F G (9, we define a positive number p(F, (9) > as 

p(F, O) := inf{|a - h{T)\^; a G \ h{0)}. 
We remark that p(-, (9) is Lipschitz continuous with Lipschitz constant 1: 

|p(F,0)-p(S,0)| <rf(F,S) (F,SgO). 
For a compact set /C C (9, we also define 

p(/C,0) :=minp(F,0). 
Let ttj = aj{T*) and 6^ = 6j(F*). Then, from the formula fl2.2p . we obtain 

= |a,„i(/i,_i(F) - + &,(/i,(F) - /i,(S)) + a,(/i,+i(F) - 

<C*d(F,S) (j = l,2,...,iV), 

where we define 

C* := max {\ai_i\ + \bi\ + \ai\}. (2.5) 

1=1,2,. ..,N 

For any F" and F^ G V* and for 6 G [0, 1], we define 

h'^ := (1 - ^)/i(F°) + ^/i(Fi) G M^. 

If there exists F^ G V* with /i(F^) = F^ is called ^-interpolation of F° and F^ The 
^-interpolation of F° G and F^ G P* is denoted by (1 - ^)F° + ^F^ := F^ G P*. 
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2.3 Polygonal motions 

We consider a moving polygon r(t) G V, where the parameter t (we call t time) belongs 
to an interval X C M. For /c G N U {0}, we call a moving polygon r(t) belongs to C^-class 



j = l,2,...,iV. 

If k > 1, we can define the normal velocity at a; G f^j{t) which is the j-th edge of T(t). 
Let nj{t) := nj{T{t)). We suppose x* G Tj{t*) and x* = {1 - e)wj^i{t*) + 0Wj{t*) for 
some 9 G (0, 1), and define x{9,t) := (1 — 9)wj_i{t) + 9wj{t) G Tj{t). Then the outward 
normal velocity of Tj{t*) at x* is defined by 

Vj{x*,t*) := x{9,t*) ■ nj{t*) = (1 - 9)Wj^i{t*) ■ nj{t*) + 9wj{t*) ■ n^it*). 

Here and hereafter, the (partial) derivative of F with respect to t is denoted by F. We 
remark that Vj{-,t) is a linear function along each Tj{t). We define the normal velocity 



where Xji'jt) ^ -^°°(r(i)) is the characteristic function of Tj{t). 

If a moving polygon r(t) belongs to a fixed equivalence class V* for all t G X, it is 
called polygonal motion in V* in this paper. Let h{t) = {hi{t), hN{t)) G be the 
height function for T{t). We remark that a polygonal motion r(t) in V* {t G X) belongs 
to C^-class if and only if h e C^il; M^), from If r(t) is a C^-class polygonal 

motion in V* , its normal velocity Vj of Tj{t) is a constant on each Tj{t) and it is given by 
Vj{t) = hj{t). We denote by VL{t) the interior domain surrounded by T{t). 

Proposition 2.2 Let T(t) he a -class polygonal motion in V* . Then we have 



on X, if the number of edges of T{t) does not change in time and Wj G C^(X; M^) for all 



of m by 



N 



V{;t) ■.= J2V,{;t)x,{;t) e L^im) 




N 



(2.6) 



Proof. From (12.21) and (12.41) . we obtain 



d 
dt 





N 
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For r G P*, the polygonal curvature of Tj is defined by 
:= — , ?7j := tany + tan^— . 

We also define the polygonal curvature of V by 

N 

The reason why this is called "curvature" is shown by the following proposition. 
Proposition 2.3 Let r{t) {t G X) he a -class polygonal motion in V* . Then 

\r{t)\ = Y,\r,{t)\^,{t)v,{t)= / ^{x,t)v{x,t)ds. 

Proof. We obtain 

N N N 

dV ^ dt 



j=i j=i j=i 

from the formula fl2.3p . I 
The polygonal curvature coincides with the crystalline curvature in the crystalline 
motion theory ([21 [11]). We, however, consider wider polygons' classes and more general 
moving boundary problems. For example, we can construct a nonconvex polygon whose all 
edges have a constant positive polygonal curvature ki = ■ ■ ■ = > (see Section [5. 2. 4p . 
We remark that such polygon is excluded in the standard crystalline theory. 



3 Initial value problem of polygonal motion 

We consider initial value problems of polygonal motions in an equivalent class. A general 
polygonal motion problem is formulated as a system of ODEs with respect to the height 
function. The notion of the constant area speed (CAS, for short) is introduced and its 
necessary and sufficient condition is given. Several concrete examples of the polygonal 
motion problems with CAS property are also presented. 

3.1 General polygonal motion problem 

We fix an equivalence class of A^-polygons V* as in Section 12.31 For an open set O G V* 
and G (0, oo], let F = (Fi, . . . , F^) be a given continuous function from O x [0, X,,) to 
with the local Lipschitz property: For arbitrary compact set /C C O and T G (0,T=k), 
there exists L(/C,T) > such that 

|F(r,t)-F(s,t)U <iv(/c,T)ci(r,s) (r,E g /c, t g [o,t]). (3.1) 
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Under the condition fl3.ip . for a compact set /C C O and T G (0,T^,), we also define 

M(/C,r) :=max{|F(r,t)|oo; T G /C, t G [0,r]} > 0. 

We consider the following initial value problem of polygonal motion. 

Problem 3.1 For a given N -polygon T* G O, find a C^-class polygonal motion T[t) G O 
{0<t<T <T^) such that 

V,{t) = F,{T{t),t) (tG[0,T], j = l,2,...,N) 

T{o) = r*. 

Under the Lipschitz condition ( 13. ip . it is clear that there exists a local solution r(t) in a 
short time interval [0,T], since Problem 13.11 can be expressed by an initial value problem 
of an ordinary differential equations for h{t). 
We consider the following assumption for Fj: 

N 

^|r,|F,(r,t) = /icAs {reO,te[o,n)), (3.2) 

where /Xcas is a fixed real number. Under the assumption (13. 2p . from the formula (12. 6p . 
any solution r(t) to Problem 13. 11 has the following property of constant area speed (CAS): 

^1^(^)1 = /"CAS- 

3.2 Examples of problems of polygonal motion 

In this section, we give some examples of polygonal motions. For several moving boundary 
problems for smooth curves, we can construct their polygonal analogues which naturally 
satisfy the basic properties such as the CAS and the curve shortening (CS, for short) 
properties. 

Problem 3.2 (polygonal curvature flow) For a given N-polygon T* G V* , find a C^- 
class family of N -polygons IJo<t<T -'^('^) ^ "P* (T < T^) satisfying 

V,(t) = -ft:,(t) (tG[0,T], J = 1,2,..., AT), 

_r(o) = r*. 

This is a polygonal analogue of the curvature flow (curve shortening problem, [3] and see 
references therein). In the theory of crystalline motion. Problem 13.21 is considered in a 
crystalline admissible class and is called the crystalline curvature motion. 

Similar to the curvature flow for smooth curves, the solution of Problem 13.21 has the 
CS property: 

.AT N 



j=l j=l 
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and the CAS property with /icAs = —2 ^^-^ tan((y9j/2): 

d ^ AT AT 

J^\m\ = -J2\^M'^^it) = -J2vj = -2 J]tan^ = const. 

j=i j=i j=i 

A numerical example will be shown in Figure [2] (left). 

Problem 3.3 (area-preserving polygonal curvature flow) For a given N -polygon 
r* G V* , find a -class family of N -polygons IJo<t<T-'^(^) C "P* (T < T*) satisfying 

V,{t) = {K{-,t)) - K,{t) {t G [0, T], J = 1, 2, . . . , iV), 

r(o) = r*. 

Here {K,{-,t)) is the mean value of k, on T{t): 

(.(, t)) = -^ [ nix, t) ds - - ' '""^^^/2) 



mUm ' ' ' \m\ \m\ 

This is a polygonal analogue of the area-preserving curvature flow (see [3] etc.). Similar 
to the area-preserving curvature flow for smooth curves, the solution of Problem 13.31 has 
the CS property: 

, V V 

d 



and the CAS property with /Icas = 0: 
d 



^^\m\ = {<;t))\m\ 



/ n{x, t) ds = 0. 

Jm 



Some numerical examples will be shown in Figure HI 

In what follows, the mean value of F on the edge Tj is denoted by 

(F),:=^ / H^)ds. 

I j\ JVj 

Let G be a bounded Lipschitz domain in M? . We define 
Og := {F G V*] n{T) D G}. 

Problem 3.4 (polygonal advected flow with constant area speed) Let us consider 
u G C^{M? \ G; R^) with div w = zn \ ^. For a given N -polygon T* G Og, find a 
G^ -class family of N -polygons IJo<t<T-'^(^) ^ (T < T^) satisfying 

'V,{t) = {u)yn, (tG[0,T], j = l,2,...,iV), 
_F(0) = F*. 



9 



The solution has CAS property with /^cas = / n ■ uds: 

JdG 

N 



/ u ■ rij ds 

,=1 -^r,{t) 



u ■ nds — I div u dx 

do Jn{t)\G 

u ■ n ds, 

do 



where n is the unit normal vector on dG pointing to interior of G. A numerical example 
will be shown in Figure [61 



4 Numerical schemes 



In this section, we propose an implicit time discretization of Crank-Nicolson type to solve 
the general initial value problem of polygonal motions (Problem I3.ip and show that it 
preserves the CAS property and has a second order accuracy. We also propose an effective 
iteration scheme to solve a nonlinear system which appears in each time step. For compar- 
isons, we also consider an explicit Euler scheme. We additionally give comments on the 
curve shortening and constant length speed properties and their numerical preservation. 



4.1 Notation 

In Section |H we consider time discretization of Problem 13.11 with the following notation. 
The discrete time steps are denoted by = to < i^i < ^2 < ■ ■ ■ < < The step size 
which may be nonuniform and their maximum size are defined by 

Tm ■= tm+1 -tm (m = 0, 1, . . . , m - 1) , T := max r^. 

0<r?i<m 

Approximate solution of T(tm) is denoted by P™ G V*. Quantities of the polygon P™ 
are denoted by ^™ = (h^^, ...,h'^):= (/ii(P™), . . . , /;.iv(r"), and := KjiT""), etc. We 
define e"^ := h{tm) - hT^ e M^. Then we have rf(P(t^),P™) = |e"^|oo. 

The discrete normal velocity = (V^", . . . , V^), which is an approximation of 
V{tm) = h{tm), is defined by 

:= (m = 0,l,...,m-l). (4.1) 

Corresponding to the formula (12.61) . the following formula holds. 



10 



This has a form of sum of areas of trapezoids and is derived from fl2.4p as follows: 



1 ^ 

|^m+l|_|^m| _ {\T']'+^\h']'+^ - \Tf\hf) 

AT 



AT ^ N 



where the last sum is equal to zero due to the equality (I2.2p . 

In the following sections, we suppose that there exists a unique solution V{t) for 
< t < T < T* to Problem 13.11 under the condition (13.11) , and that discrete time steps 
Q = tf) < ti < t2 < ■ ■ ■ < tfh < T are given a priori such as the uniform time stepping 
tm = 1TLT. We adopt the uniform time increment in the numerical examples in Section[5l It 
is, however, possible to apply any a posteriori adaptive time step control scheme. Similar 
to the finite time extinction of the curvature flow of smooth curves, even in the polygonal 
motions, the solution polygon often has singularities in finite time. For instance, |rj(t)| 
tends to zero, in other words, tends to infinity. A posteriori adaptive time step 

control will be required near the blow-up time for accurate computation. 

4.2 Second order implicit scheme 

We consider the following implicit scheme for Problem 13. II 

Problem 4.1 For a given N -polygon G (9 and given time steps = to < < ^2 < 
■ ■ ■ < tm ^ T , find polygons F™ E O {m = 1,2, ... ,Th) such that 

yr = (r™+'/'> W1/2) (m = 0, 1, 2, . . . , m - 1, J = 1, 2, . . . , AT), 

r° = F„ 

where r™+-'^/^ and tm+1/2 are the 1/2-interpolations: 



pm _|_ pm+l T- 



This is a generalized version of the scheme presented in [13] for area-preserving crystalline 
curvature flow. 

Theorem 4.2 We suppose the CAS property (13.21) . Let F™ G O (m = 1,2, . . . ,m) be a 



solution of Problem 4^.1. Then it satisfies 

|^m+i| ^ (m = 0,l,...,m-l). 

In other words, \VL'^\ = \Q{tm)\ holds if the exact solution Q(t) of Problem \3.1\ exists. 
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Proof. Since \Tf+^^^\ = i\Tf\ + \Tf+^\)/2, we have 



E|pm+l/2|p/pm+l/2 . , ^ _ ,, 

|ij l-Tjli , rm+l/2j — /^CAS, 



Tm ■ -, 

from the formula (14. 2p . I 
We remark that the numerical scheme Problem 14.11 inherits the CAS property but 

does not depends on the area speed /Xcas- 

Since Problem 14.11 is an implicit scheme, it is not clear whether r™+^ G O can be 

determined uniquely from the previous polygon P™ G O, the time tm, and the time step 

size Tm- Another question is how to solve the equations 

/pm _|_ pm+1 \ 



to obtain (approximation of) P"^+^ numerically. The answers to these questions will be 
given in Theorem 14.41 and 14. 5[ 

We fix P G O and i G [0, T), which correspond to P™" and tm+i/2, respectively. Let /C be 
a compact convex set in V* with f G /C C C For S G /C and f G (0, p(f, C)M(/C, T)-i), 
we can define S G C by 

h{t) = h{t) + fF {^^J 

We define A(S) := A(S; f , t, f) := S. Then A becomes a mapping from K, to (9. We have 
the following lemma. 



Lemma 4.3 Let e G (0,p(f,(9)) and A G (0,1) he fixed, and let fC := B{V,£). Suppose 
that T satisfies the condition: 

,^ ' ^ 2A 
< f < min <^ T - t, 



M(/C,T) L(/C,T) 
r/ien A ma]9S /C mto /C and satisfies 

d{X{T}),K{T?)) < \d{T}, T?) (S\ G /C). (4.4) 
Namely, A is a contraction mapping on /C and there exists a unique fixed point of A inKL. 
Proof. Let S G /C. Since /C is convex, (f + S)/2 G /C holds. We also have 



h{t) - h{t) 



< fM(/C,T) < e. 
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This estimate shows that A is a mapping from /C into itself. The estimate fl4.4l) is proved 
as follows: 



d(A(E'),A(E=)) 

= |MA(E'))-MA(E"))|oo = f 



< TL{t,T)d 
= fL(/C,T) 



/f + f + S2 
^(f) + ^(Si) /i(f) + h(E2) 



2 2 

L{t,T)d{T},Y?) < \d{T},Y?) . 



The following theorem gives us an efficient numerical scheme to obtain r™"*"^. The proof 
is clear from Lemma [4.31 

Theorem 4.4 Let K he a compact set in O and let e G (0,p(/C, O)). We define 



/C, := |J^(^'^)- 



se/c 



For fixed m (< m) m Problem \4-l\ we assume that T"^ G JC and 

e 2A 



Tm < min 



M(/C„T)'L(/C„T) 



where A G (0, 1). T/ien t/iere exists uniquely T"^'^^ g i?(r™,£:) satisfying (14.31) . 

Furthermore, T"^'^'^ is a fixed point of the contraction K„i '■= A.(- ; F™, t^, x^) in BiV™-, e), 
and is given by the limit o/ A,^(r'") as u —>■ oo with the following estimate 



d(r™+\ Aj;(r'")) < A''d(r™+\ r™) (z/ g n). 

From this theorem, A'^{r'^) for sufficiently large u gives a satisfactory approximation 
of F'"^^. An iteration algorithm based on this idea will be given in Section 15. 1[ The 
following theorem shows the second order convergence of the implicit numerical scheme 
(Problem iJ]). 

Theorem 4.5 We suppose that {r(t)}o<t<T be a C^^^-class solution of Problem \3. 1\ for 
= 0, 1, or 2. There exists 5* > 0, r* > 0, C > and a non- decreasing function uj{a) > 
with 



ujia) 



' o{a^) if k = or I, 



as a I 0, 



(4.5) 
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such that, z/ci(r*,r°) < 5* and r < t* , then T"^ e O {m 



1,2, 



m) are inductively 



determined by the implicit scheme (Problem \4. 1^ and 



max d{T{tm), T"^) < cj(r) + Cd{T{0), r°), 

0<m<m 

holds. 

Proof. We put p := p({r(t); < t < T},0), and fix 5 G (0,p) and e E (0, p - 5). We 
define 



/C:= U 5(r(t),5), /C,:= |J^(S'^)' 

o<i<r se/c 

(7- \ — 1/a 
(0<a<2/L), 



Pr. 



+ 



[m = 0, 1,2, . . . ,m). 



where a non- decreasing function Ci;(a) (0 < a < T), which satisfies (14. Sp . will be defined 
in (14 .Qp later. Since R{-) is an increasing function, there exists 6* > and r* > such 
that 



i?(r 



e 



For m = 0, 1, 2, . . . , m — 1, we will prove the following inductive conditions: 

r"^G/C, Pm<R{ry-po ^r^+^G/C, < i?(r)*-+>o. (4.6) 

The conditions F*^ G /C and po ^ Ri'^)^Po for the case m = are obviously satisfied. 

Let us assume the conditions F'" G /C and pm < Ri^Y'^Po for a fixed m. Then, from 
Theorem 14. 4[ there exists F™"''^ uniquely in B{T"^,e) C /C^, and we have 

^m+l _^m^ _ _ ^^ym ^ ^^m ^ {y{t^^,,^) - F™) } , (4.7) 

j-m h{tm+i) — h{t„i 



- hit. 



m+l/2) 



The last term of (14. 7p is estimated as follows. Since F(tm+i/2) E K, C JCe and r™+^/^ G 



5(F'^,£:) C /C^, we have 

|V(Wl/2)-V^'"l 



|F(F(t, 



^m+1/2) — F(F'""^^''^, t. 



m+l/2j 



<Ld{T{t^+u2),r^^'^') 



h{t. 



m+l/2 ) 



(4. 



where 



.m _ Htm) + fa(Wl) , ^ 
S. •— "'lrm+l/2j 
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Combining fl4.7p and fl4.8l) . we obtain 



loo 



< w 



< 1/3™ I -I- T I /J"* I -I- T- J 

I '^m^i /I m+1 



-(e'" + e'"+^) -C 



e-+i|^ + |e-U)+r^(|r|oo + i:|C 



By the Taylor expansion, we can obtain an non-decreasing function co'(a) (0 < a < T) 
which satisfies the condition (14.51) and the inequahty 



irioo+i:inoo<cu(r). 

Hence, we have 

Trnlj 



TmL 



and this inequahty is equivalent to 



oo + TmUjiT) 



Tm.L 



Prn+l < 1 + 



TmL 



Pr, 



From the inequalities 



1 



> 1 

2 / - V 2 



t,„/t 



and {l + ^]< e^-^/^ 



we obtain 

Pm+l < I 1 



The condition r'"^^ G /C follows from this estimate as 



^m+l I 



a; 



Hence, we have proved (14.61) . which leads us to the estimate: 



*\T / 1^0 



*\T 



(4.9) 



eloo + ^ ) - ^ ^ i?(r*)^|e°U + ^^^^^^ M^)- 



The assertion of the theorem is obtained by putting C := R{t*)^ and denoting the last 
term L~^{R{t*Y — again by uj{t). I 
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4.3 Euler scheme 



For Problem I3.H one of the simplest numerical scheme is the following explicit Euler 
scheme: 

Problem 4.6 For a given N -polygon G (9 and time steps = to < < ^2 < ■ ■ ■ < 
tfh < T , find polygons F™ E O (m = 1, 2, . . . , m) such that 

Vp = F,{T"^,t^) (m = 0,l,2,...,m-l, j = 1, 2, . . . , iV), 

ro = F,. 

The explicit Euler scheme is simple but it has only first order accuracy. In particular, 
for polygonal motions with CAS property, we are required to use a more accurate scheme 
such as Problem 14.11 in order to keep its CAS property numerically. Similarly to the case 
of the implicit scheme (Theorem 14. 5p . the convergence theorem of the Euler scheme is 
stated as follows. 

Theorem 4.7 We suppose the condition (13. ip and that {F(t)}o<t<T be a C^^^ -class so- 
lution of Problem \3.1\ for k = or 1. There exists 6*>0, t*>0,C>0 and a 
non- decreasing function uj{a) > with 

fo(l) ifk = 0, 

uj{a) = < as a 10, (4.10) 

[0{a) tfk = l, 

such that, if d{T* , F") < 5* and r < r* , then F™ G O (m = 1, 2, . . . , m) is determined by 
the Euler scheme (Problem \4.0^ and satisfies the estimate 

max d(F(t„), F"') < cu(r) + Cd(F(0), F°). 

0<m<m 

Proof. We define 

r := ^^'"^^^ ~ ^^'"^ - ktr^) (m = 0,...,m-l). 

Then, by the Taylor expansion, we are able to find an non- decreasing function uj{a) 
(0 < a < T) which satisfies the condition (I4.10p and the inequality 

ir loo < uj{r). (4.11) 

We put p:= p ({F(t); < t < T}, C), and fix 5 G (0, p) and e G (0, p - 5). We define 



/C:= U i?(F(t),5), /C,:= \jB(j:,e), 

o<t<T se/c 

L:=L(/C,T), p^:=|e'"U + ^ (m = 0, 1, 2, . . . , m). 
There exists 6* > and r* > such that 



M(/c,r) 
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For m = 0, 1, 2, . . . , m — 1, we will prove the following inductive conditions: 

r'^G/C, p^<e*'"% ^ ^r^+^G/C, < e*-+^^po. (4.12) 

The condition for m = are obviously satisfied. 

Let us assume the conditions F"^ G /C and pm < e*™^po for a fixed m. Then, from the 
condition r* < e/M{}C,T), F^+i belongs to B(T"',e) C /C^, and we have 

Since F(tm) and F™ both belong to /C, we have 

- ^'"loo = \F{T{tJ,tJ - ^(F^OI^ < Lrf(F(U,F'") = L|e"^|^.(4.14) 
Combining (HATII . (H13|) and (SHI), we obtain 

|e'"+^|oo < (1 + r^L) |e"^U + rmOo{r), 

and 

Pm+i < (1 + rmL)pm < e"'"^(e*'"^ po) = e^^+^^po. 

Since 

the condition F™+^ G follows. Hence, we have proved (14.121) . which leads us to the 
estimate: 

|e"U < e'-- (|e«U + ^) - ^ < e-|eV + ^-(-)- 

The assertion of the theorem is obtained by putting C := e^^ and denoting the last term 
2^-1 ^gLr _ i)(^(t) again by L<j(r). I 

4.4 Curve shortening and constant length speed property 

As seen in Section 13.21 many moving boundary problems hold the CS property: 

>)l<o. 

From (12. 3p . a necessary and sufficient condition for CS property is 

N 

J2v,F,{T,t) < (F G O, t G [0,T,)). (4.15) 
i=i 

Similarly to the CAS property (13.21) . we can also consider constant length speed (CLS, 
for short) property: 

^|F(t)| = /icLS. 
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A necessary and sufficient condition for CLS property is 

N 

Y,V,F,{T,t) = /icLs iTeO,te [0,T,)). (4.16) 
An example with CLS property is the constant speed motion: 



N 



F^{T,t) = l (j = l,...,Ar), /ieLs = 2^tan2 



F,(r,t) = f^=^'/' L -^^T) {j = l,...,N), /icLs = 0. (4.17) 



Another example is the length-preserving polygonal curvature flow: 

Ef=iir,ls-(r)^ 

It is easy to check that both the second order implicit scheme (Problem 14. ip and the 
explicit Euler scheme (Problem 14. 3p inherit the CS and CLS properties. Namely, under 
the condition (14.151) . we have 

|pm+l| < |pm| (m = 0,l,...,m-l), 

and, under the condition (I4.16p . we have 

|r-+i| = IF'"! + /icLsT^ (m = 0, 1, . . . , m - 1). 

A numerical simulation for the length-preserving polygonal curvature flow will be shown 
in Figure [Hi 



5 Numerical computation 

We describe an algorithm of our second order implicit scheme and show some numerical 
results. In this section, F™ G V* {m = 0,1, ... ,m) denotes the numerical solution 
computed by the algorithm described in Section 15.11 All computations are performed in 
double precision. 



5.1 Algorithm 

We describe a numerical procedure of Problem 14. 1[ We suppose that an initial N- 
polygon F° = IJ^i is given in a prescribed equivalence class V*, i.e., V* = P[F°]. 
In other words, the set of normal vectors {rij}^^ for V* and the set of heights of F° 
/i° = (hi, . . . , /i^) G are given. The outer angles {ipj}jLi and the quantities {aj}jLi, 
{bj}^^i are computed from from {n^}^^. We fix the maximum computation time T^, and 
the uniform time step r = T*/m with the maximum time step m. 

Then F™+^ G V* is determined successively from F*" G V* at the m-th discrete 
time tm = 1TLT for m = 0, 1, ... ,m — 1 as follows. We suppose the set of heights of 
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pm j^m ^ ^^m^ j^m^ ^ -^N given. We can calculate the j-th vertex wj^ of 
by (12.11) (j = 1,2, ...,A^). Our algorithm including the iteration scheme to obtain an 
approximation of r"^~^^ is as follows. 

(1) Put h := /i™. 

(2) Define t e V* with ^(f) = h and put h := h. 

(3) Compute h := h"" + F(f,t„+i/2) r/2. 

(4) If \h — h\oo < s/2, then go to step (6). 

(5) Go to step (2). 

(6) Put /i^+i := 2h - h"". 

We note that f and T (with h{T) = h) in step (3) correspond to (A^(r'") + r'")/2 
and (A5:^^(r™) + r™')/2, respectively. In the stopping condition (4), we choose a small 
parameter e > 0. In the following numerical computations, we took e = 10^^^. 

5.2 Numerical examples 

In the following examples, several numerical computations of the evolution of A^-sided 
polygons will be shown. The numerical solutions were computed until the time with 
the uniform time increment r = T^/rh, where m is the maximum time step. The figures 
are depicted every M-th time step. The problems except Example 7 have the CAS 
property with /Icas; and the numerical solution keeps this property with the error A = 
maxo<m<m lAtcAs ~/^cAsl; where /i^g = (|f2'""'"^| — |f2™|)/r is the m-th discrete area speed. 
The problem in Example 7 has CLS property with /Xcls, and the numerical solution keeps 
this property with the error A = maxo<m,<m |/^cls ~/^ct,sI; where /i^g = (IP^+^I — iF^D/r 
is the m-th discrete length speed. The following two tables indicate the data A^, T*, r, M 
and A in each example. 





Ex.1: Figure E] 


Ex.2: Fi 


gure[3] 




(left) 


(middle) 


(right) 


(left) 


(right) 


N 


5 


7 


22 


7 




0.2801 


0.3136 


0.335 


1.55 


1.55 


T 


10-6 


10-4 


10-7 


M 


2801 


3136 


3350 


775 


775000 


A 


5.04 X 10-^° 


7.62 X 10-^° 


1.56 X 10-^5 


2.87 X 10-^^ 


2.96 X 10^ 



Table 1: Numerical parameters and A for Examples 1 and 2. 
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Ex.3: Figure H 


Ex.5: Figure [6] 


Ex.6: Figure [7| 


Ex.7: Figure El 






(upper) 


(lower) 


(upper) 


(lower) 




N 


9 


12 


12 


32 


18 


% 


7.56 


19.4 


20 


10 


3.65 


0.27 


T 


10 


-5 


10-4 


10-4 


10-4 


M 


37800 


97000 


10000 


5000 


1825 


27 


A 


1.51 X 10-9 


1.07 X 10-9 


2.81 X lO-'^ 


4.26 X 10-9 


1.42 X 10-10 


5.11 X 10-11 



Table 2: Numerical parameters and A for Examples 4-7. 
5.2.1 Example 1 — polygonal curvature flow 

Figure [2] indicates the evolution of solution polygons to Problem 13.21 starting from the 
initial polygon being the outermost iV-sided polygon which is a combination of an upper 
half of a regular 2(iV — 2)-polygon and a triangle. Each solution polygon evolves from out- 
side to inside and has the CAS property with /Icas = ~2 ^^-^ tan((y9j/2). The numerical 
solutions keep the CAS property very accurately as shown in Table [H 




Figure 2: Evolution by polygonal curvature flow. 



5.2.2 Example 2 — backward polygonal curvature flow 

Problem 13. 21 can be computed backward in time. Figured (left) indicates the evolution of 
solution polygons to the backward polygonal curvature flow Vj{t) = Kj(t) (j = 1, 2, . . . , 7). 
The initial polygon is the innermost 7-sided polygon and the solution polygons evolve 
from inside to outside. The above process can be followed by our second order scheme 
accurately. 

We note that the backward curvature flow for smooth curves is ill-posed since it 
becomes a backward parabolic problem. Actually, even in the case of 7-sided polygon's 
motion, it is hard to compute the backward polygonal curvature flow by using the Euler 
scheme (Problem I4.6p . Figure [3] (right) indicates an easy breakdown of the Euler scheme 
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in spite of using a smaller r than the one in the second order scheme. 




Figure 3: Simulations of the backward polygonal curvature flow by the second order scheme 
(left) and by the Euler method (right). 



5.2.3 Example 3 — area-preserving polygonal curvature flow 

Figure H] (middle-upper/lower) shows two examples of polygonal motions according to 
Problem 13.31 The initial polygons are given as Figure H] (left-upper/lower). Figure H] 
(right-upper/lower) shows the final polygon and the initial polygon (dotted curve). The 
solution has CAS property with /Icas = 0. 

In both upper and lower examples, there exist stationary solutions as shown in Figure [5] 
(middle/right). The polygon starting from a symmetric initial shape approaches to one 
of the stationary solutions and stays there for a while. However, since the stationary 
solution has a saddle-point instability, after a while, the polygon is drifted away from the 
stationary solution along the unstable manifold and loses its symmetry. 




Figure 4: Evolutions by the area-preserving polygonal curvature flow. 
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5.2.4 Example 4 — 



stationary solutions 



A polygon which has a constant polygonal curvature (i.e. Ki = ■ ■ ■ = Kn) is a stationary 
solution of Problem l3.3[ Obviously, regular polygons are stationary solutions. Besides the 
regular polygons, we have infinite many stationary solutions. For instance, an n-fold star 
shaped polygon is a stationary solution as well as the 6-fold star (Figure O (left)), and an 
n-fold non-sharp star shaped polygon is also a stationary solution as well as the 3-fold/4- 
fold non-sharp stars in Figure [5] (middle/right). For the ra-fold non-sharp star polygon, 
there are two kinds of outer angles ipQ < and ipi > and two kinds of edge lengths 
di < d2 with the corresponding polygonal curvatures: Ki = (tan(y9o/2) + ta.n{(pi/2))/di 
and ^2 = 2tan{(pi/2)/d2. We have constant polygonal curvature polygon if di/d2 = 
(tan(y9o/2) + tan(y9i/2))/(2tan((y9i/2)). The polygon in Figure [5] (middle/right) belongs 
the same equivalence class of Figure H] (upper/lower). 



5.2.5 Example 5 — polygonal advected flow with constant area speed 

Figure [6] (middle) shows an example of polygonal motions according to Problem 13.41 with 
u{x) = x/ (27t\x\'^) which is a divergence-free vector field defined on \ {0}. The initial 
polygons are given as Figure E] (left) whose center is the origin. Figure E] (right) shows the 
final polygon and the initial polygon (dotted curve). The problem has the CAS property 
with /icAs = 1- The numerical solution keeps the CAS property very accurately as shown 
in Table H 




Figure 5: Stationary solutions of Problem 13.31 




Figure 6: Evolution by polygonal advected flow with constant area speed. 
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5.2.6 Example 6 — area-preserving polygonal advected-curvature flow 

Figure [7] (middle-upper/lower) sliows two examples of polygonal motions according to 
combination of Problem 13.31 and Problem 13.41 i.e., 

= t)) - ^jit) + ■ {j = l,...,N). 

The divergence-free vector field is given by u{x) = xiX2(— xi, a;2) (upper) and u{x) = 
{—Xi,X2) (lower), respectively. The common initial polygon is given as Figure [7] (left- 
upper/lower), where the center is the origin the vertices are on the ellipse with ratio 
3:1. Figure [7] (right-upper/lower) shows the final polygon and the initial polygon (dotted 
curve). The both problems have CAS property with fi^As = (area preserving), and the 
numerical solutions preserve their areas very accurately as shown in Table [2j 





Figure 7: Evolution by area-preserving polygonal advected-curvature flow. 



5.2.7 Example 7 — length-preserving polygonal curvature flow 

Figure E] (middle) shows the evolution of solution polygons to length-preserving polygonal 
curvature flow fl4.17p given in Section 14.41 The initial polygon is the 18-sided polygon 
in Figure [8] (left). Figure [8] (right) indicates the final polygon and the initial polygon 
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(dotted curve). At the time close to T^,, the length of an edge (pointed by the arrow) 
tends to zero, and the computation stops. The length-preserving property (/icLs = 0) is 
numerically realized with high accuracy as shown in Table [2j 




Figure 8: Evolution by the length-preserving polygonal curvature flow. 
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